Finite-size effects in BaTiOs nanowires 
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The size dependence of the ferroelectric properties of BaTi03 nanowires is studied from first- 
principles. We show that the ferroelectric distortion along the wire axis disappears below a critical 
diameter of about 1.2 nm. This disappearance is related to a global contraction of the unit cell 
resulting from low atomic coordinations at the wire surface. It is shown that a ferroelectric distortion 
can be recovered under appropriate tensile strain conditions. 



Ferroelectric oxides constitute an important class 
of multifunctional compounds, attractive for numerous 
technological applications*- Pushed by the miniatur- 
ization constraints imposed by future devices, there is 
now an increasing interest in the synthesis and use 
of nanosized ferroelectrics 2 ". Nowadays, various ferro- 
electric nanostructures - ultrathin films*^, nanowires 5 , 
nanotubes^i, nanoislands 8 - can be grown with control 
at the atomic scale. However, bringing ferroelectrics to 
the nanoscale raises fundamental questions. 

Ferroelectricity is a collective effect, resulting from 
a delicate balance between short-range and long-range 
Coulomb interactions. In confined structures, both in- 
teractions are modified with respect to the bulk and it 
is commonly believed that ferroelectricity is altered and 
eventually totally suppressed when the system reaches a 
critical size sometimes referred to as the "ferroelectric cor- 
relation volume'—. First-principles simulations are a fun- 
damental tool to understand size effects in ferroelectrics 9 . 
Recent first-principles studies provided valuable insight 
into the behavior of thin films and highlighted the crucial 
role of the depolarizing fieldifi.ilii£. A few works based on 
effective Hamiltonians»li««*ii and shell models!*- also con- 
sidered other low dimensional systems, such as nanorods 
and nanodisks. However, despite these efforts, very little 
is known on the evolution of ferroelectric properties in 
nanowires with respect to their size and shape. 

In this letter, we investigate from first-principles the in- 
fluence of finite-size effects on the ferroelectric properties 
of BaTiC-3 nanowires. We identify a critical diameter be- 
low which the ferroelectric distortion along the wire axis 
is totally suppressed. We assign this behavior to surface 
atomic relaxation effects. Below the critical diameter, 
a ferroelectric state can be recovered under appropriate 
tensile strain conditions. The crucial role of the longi- 
tudinal interatomic coupling along Ti-0 chains and the 
weak interchain interactions are also discussed. 

We focus on stoichiometric infinite BaTiC-3 nanowires. 
As illustrated in Fig. EJ we consider model systems built 
from individual BaTiC>3 clusters that are assembled into 
infinite chains, themselves combined into wires of increas- 
ing diameters. These stoichiometric wires have two BaO 
and two TiC>2 lateral surfaces, as displayed in Fig. Eft). 
This asymmetry is responsible for a dipole perpendicular 
to the wire axis. However, it will not be further discussed 
since it is not relevant to the present study^. In what fol- 




FIG. 1: (Color online) BaTi03 stoichiometric wires: BaTiOs 
clusters (a) assemble into infinite chains (b), themselves gath- 
ered into nanowires (c) . The diameter of the nanowires is de- 
fined from the number of Ti-0 chains (n) assembled together 
within the wire. Panels (c) and (d) correspond to n = 4. 



lows, we will focus on the eventual disappearance of the 
ferroelectric instability along the wire direction z. The 
equilibrium lattice parameter along z will be referred to 
as ao, and the wire diameter will be specified by an inte- 
ger n that corresponds to the number of Ti-0 chains that 
are gathered together within the wire. The alternating 
Ti-0 distances along z will be called d\ and d,2, as shown 
in Fig. Eft). 

Our calculations have been performed within the lo- 
cal density approximation (LDA) to density functional 
theory (DFT) using a numerical atomic orbital method 
as implemented in the Siesta codeii. Technical details 
are similar to those of Ref. Some of the calculations 
have been double checked using the ABINIT packageii". 
Periodic replicas of the wires are separated by more than 
10 A of vacuum. We explicitly checked that the replica 
of the wire generated by periodic boundary conditions 
is not affecting our results 2 ^. The atomic positions have 
been relaxed until the maximum component of the force 
on any atom is smaller than 0.01 eV/A. 

For the thinnest wire (n = 1, Fig. Eft)) the bond 
lengths and the equilibrium lattice constant ao of the 
relaxed structure reported in Table [J are significantly 
shorter than in bulk. This can be explained from the 
low coordination of the different atoms in the structure, 
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TABLE I: Interatomic distances and equilibrium lattice pa- 
rameter a along the wire axis, for the n = 1 wire. 0/ (resp. 
On) refers to oxygen atoms along the Ti-0 (resp. Ba-O) 
chains. Results for the bulk cubic paraelectric unit cell of 
BaTi03 are shown for comparison. Units in A. 





Ti-Oj 


Ba-Oi/ 


Ba-Ti 


ao 


Wire 


1.88 


2.52 


3.07 


3.60 


Bulk 


1.97 


2.78 


3.41 


3.94 



the lack of some interactions being compensated by the 
strengthening of the the remaining bonds. The system 
remains insulating with a bandgap of 3.3 eV. The relaxed 
structure presents a plane of symmetry perpendicular to 
z (di — da) and is therefore not spontaneously polar- 
ized along the wire direction, demonstrating that ferro- 
electricity has disappeared. However, imposing a tensile 
strain, we can recover a polar structural distortion along 
the Ti-0 chains (d\ ^ cfe) above a critical lattice con- 
stant a c = 4.05 A. This strain induced phase transition 
is second order. Above a c , both di — ^2 and the polar- 
ization, calculated from the Berry phase approach jSli 2 ^ 
increase linearly with ao- 

The fully relaxed n = 4 wire, shown in Fig. QJc), is 
also non-polar along z. However, a c has been reduced to 
3.88 A and ao is now larger (3.80 A) as expected from 
the increase of the mean atomic coordination numbers 
with respect to the n — 1 wire. Figure El represents the 
evolution of ao and a c for wires of increasing n. We ob- 
serve that the difference, a c — ao, progressively decreases. 
A crossover between both curves is observed at n c = 9, 
that corresponds to a diameter of 1.2 nm. Above n c , 
unconstrained wires are ferroelectric whereas below this 
value, a spontaneous polarization only appears under ten- 
sile strain. 
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FIG. 2: Evolution of a c (white diamonds) and ao (black di- 
amonds) as a function of the number n of Ti-0 chains in a 
transverse section of the wire. 

In Fig. [2j a c keeps a nearly constant value around 
3.85 A (except for n — 1) while the disappearance of 
the ferroelectric distortion at small sizes originates from 



the progressive decrease of ao. The ferroelectric behav- 
ior of nanowires appears therefore to be monitored by 
the strong sensitivity of ferroelectricity to the unit cell 
volume, the latter affected by low coordination effects at 
free surfaces. As n increases, we recover a polar distor- 
tion, however, its amplitude only progressively increases 
as a slowly evolves to its bulk value. 

Similarly, below n c , a second-order phase transition to 
a polar state is observed under tensile strain, but the 
resulting polarization strongly depends on the lattice pa- 
rameter. This is illustrated in Fig. for n — 4, where 
we quantify the polar distortion through the anisotropy 
of the interatomic distances (d2 — di) along the distinct 
Ti-0 chains. This provides a qualitative measure of the 
contributions to the polarization of the individual chains 
since the Born effective charges (along z) are reasonably 
similar from one chain to the other (variations smaller 
than 20 The strong inhomogeneity of the distor- 

tions for the different Ti-0 chains highlights a deep in- 
fluence of the surface: the ferroelectric distortion paral- 
lel to the wire surface is enhanced at Ti02 surfaces and 
reduced around the BaO surfaces, a feature similar to 
what was previously reported for free-standing slabs^. 
Note that the transition to the polar state does not oc- 
cur at the same lattice constant in the different chains, 
giving rise, within a given range of lattice constants, to 
mixed ferroelectric-paraelectric states in these nanosized 
systems. 
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FIG. 3: Local distortions as a function of the lattice parameter 
for the n = 4 wire. Labels of the chains as in Fig.^d). 

The distinct behavior of the different Ti-0 chains 
within the same wire suggests weak inter-chain interac- 
tions. This is further confirmed by the fact that the po- 
larization of each chain can be reversed almost indepen- 
dently, with very small changes in the total energy of the 
system and in the amplitude of the polar displacements. 
It is the longitudinal coupling of the atomic displace- 
ments inside individual Ti-0 chains that appears as the 
main origin of the ferroelectric distortion. This behavior 
is similar to what was previously reported for the bulk 2t> . 
In the phonon dispersion curves of BaTiO^, the unsta- 
ble ferroelectric mode at T remains similarly unstable at 
X and M but not at R, pointing out a strong anisotropy 
of the interatomic force constants 25 , a significant cou- 
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pling of the atomic displacement along Ti-0 chains, and 
the absence of strong interchain interactions. The result- 
ing "chain- like" character of the ferroelectric instability 
is intrinsic to BaTiC>3. It is preserved in nanowires and 
provides a general argument to understand the survival 
of ferroelectricity along the wire axis down to very small 
radii. The situation might be different in other com- 
pounds such as PbTiC>3 in which the ferroelectric insta- 
bility is more isotropic^. 

In summary, the ferroelectric behavior of infinite sto- 
ichiometric BaTiC>3 nanowires has been studied from 
first-principles. In BaTiC>3, the ferroelectric instabil- 
ity exhibits a marked "chain- like" character so that, in 
nanowires, it could a priori be preserved down to very 
small sizes. Nevertheless, at small radii, low atomic co- 
ordinations at the surface produce a contraction of the 



unit cell that is responsible for the suppression of the fer- 
roelectric distortion at a critical radius estimated around 
1.2 nm. Below this radius the ferroelectric distortion has 
disappeared at the equilibrium volume but can be recov- 
ered under appropriate tensile strain conditions. 
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